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Abstract 

Casimir effect in the planar setting is described using the boundary state for- 
malism, for general partially reflecting boundaries. It is expressed in terms of the 
low-energy degrees of freedom, which provides a large distance expansion valid for 
general interacting field theories provided there is a non- vanishing mass gap. The 
expansion is written in terms of the scattering amplitudes, and needs no ultraviolet 
renormalization. We also discuss the case when the quantum field has a nontrivial 
vacuum configuration. 

1 Introduction 

The Casimir effect can be considered as a response of the ground state in a quantum field 
theory to the presence of boundary conditions. Therefore it is natural to seek a relation to 
the approach known as boundary quantum field theory started in two-dimensional space- 
time by the seminal paper of Ghoshal and Zamolodchikov [T]. Recently we have developed 
and extended this formalism to quantum field theories in arbitrary space-time dimensions 
and applied it to the Casimir effect 02, [H [5] . Here we give a short review of our results. 



2 Boundary state formalism 

2.1 The concept of the boundary state 

Following [5] we consider an Euclidean quantum field theory of a scalar field $ defined in 
afl + 1 dimensional half space-time, parameterized as (x < 0,y,r), in the presence of a 



1 



codimension one flat boundary located at x = 0. The correlation functions defined as 



y u n) . . . $(xjv, yjv, = j p$ e _ s[$] 

contain all information about the theory. The measure in the functional integral is provided 
by the classical action 

S[$] = fdrf°°dy\ [ dx f^(V$) 2 + [/($)) + f/ B ($(x = 0,y,f)) 

which determines also the boundary condition via the boundary potential Ub- 



d $1 — 



2=0 



Here we assume for simplicity that the boundary term does not depend on the temporal 
(i.e. y) derivative of $, which means that there are no boundary degrees of freedom with a 
temporal dynamics independent of the bulk (it may depend on derivatives with respect to r, 
which is the reason for the variational derivative 5). The bulk interaction U is constrained 
by the requirement that the bulk spectrum must possess a mass gap m. 

This Euclidean quantum field theory can be considered as the imaginary time ver- 
sion of two different Minkowskian quantum field theories. We can consider t = —iy as 
Minkowskian time and so the boundary is located in space providing nontrivial bound- 
ary condition for the field <3>. The space of states in this Hamiltonian description is the 
boundary Hilbert space TCb determined by the configurations on the equal time slices. Hb 
contains multi-particle states and is built over the vacuum state, obtained in the presence 
of the boundary condition (|0)#), by the successive application of particle creation oper- 
ator^. In the asymptotic past the particles do not interact and behave as free particles 
travelling towards the boundary; thus 

Hb = {afnih^h) . . .af n (k Nj k N )\0} B , h > . . . > k N > o} 

where the operator af n (k, k) creates an asymptotic particle of mass m with transverse (i.e. 
x-directional) momentum k and parallel (i.e. parallel to the boundary) momentum k. 



The corresponding energy is u(k, k) = y k 2 + k 2 + m 2 = Jk 2 + m e ff(k) 2 where m e s(k) 



k 2 + m 2 is the effective mass of a particle with parallel momentum k as seen in the two- 
dimensional space-time formed by t and x. Instead of k, we shall also frequently use the 
rapidity parameter •& defined by 

to = m e s(k) coshi? , k = m e s(k) smh.-& (2.1) 



1 One can also introduce particle-like excitations confined to the boundary [3] ('surface plasmons'), but 
for simplicity we do not consider them here. 
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Figure 2.1: The two Hamiltonian descriptions, with a representation of the amplitudes R 
and K 2 . 



In the Heisenberg picture the time evolution of the field 



$(x,t,r) = e iHBt $(x,0,r)e 



-iH B t 



is generated by the following boundary Hamiltonian 



I 



df 




dx Qn t 2 + l -{d x §Y + 1 -0§Y + [/($)) + U B (*(x = 0))] (2.2) 



The correlator can then be understood as the matrix element 

($(xi,yi,fi) . . . ${x N ,y N ,r N )) = B (0\T t ($(xi,ti,fi) . . . §(x N , t N , f N )) \0) B 

where T t denotes time ordering with respect to time t, and the vacuum |0)_b is normalized 
to 1. 

The formulation of asymptotic states and fields, together with the relevant reduction 
formulae (which generalize the LSZ approach to boundary QFT) was given in [21 E]. In 
[3j we also gave the appropriate generalization of Landau equations, Coleman-Norton in- 
terpretation and Cutkosky rules, together with an example of one-loop renormalization of 
boundary interaction (where we considered the case of sine-Gordon model in two space- 
time dimensions). Elastic reflection of a particle from the boundary (see figure |2~T1) is a 
process with one particle of energy ui and parallel momentum k both in the incoming and 
outgoing stated, whose transverse momentum k changes sign. Its amplitude is the reflec- 
tion factor R(u,k) which can only depend on \k\, as a result of rotational invariance in 
the directions parallel to the boundary; it is not necessarily unitary due to the possible 
existence of inelastic processes. 

Alternatively we can also consider r = — ix as Minkowskian time, as depicted in figure 
12.11 In this case the boundary is located in time and we can use the usual infinite volume 
Hamiltonian description. The Hilbert space is the bulk Hilbert space H spanned by multi- 
particle in states 



2 Energy and parallel momentum are conserved due to the unbroken translation invariance in the direc- 
tions parallel to the boundary. 



H = {Af(K h h) . . .A+(K N ,k N )\0) 



kl > . . . > fcivj 



where k is the momentum in the y direction, and the energy corresponding to the time di- 
rection is given by lu(k, k) = \J m 2 + k 2 + k 2 . One can again use a rapidity parametrization 
in this channel defined by 

k = m e g(k) sinh'i? , u = m e g{k) cosh^ (2.3) 

Time evolution 

$(T,y,r) = e iHT my,r)e- iHT 
is generated by the bulk Hamiltonian 

H = J drj^dy + ^(c\$) 2 + \{d^f + [/(<&)) (2.4) 

and the boundary appears in time as a final state in calculating the correlator: 

y u n) . . . <5>(x N , y N , r N )) = (B\T T ($(n, y u n) . . . $(t n , y N , f N )) |0) 

The state (B\ is called the boundary state, which is an element of the bulk Hilbert space 
and is defined by the equality of the two alternative Hamiltonian descriptions 

(B\T T ($(n, y u n) . . . $(tjv, y N , r N )) |0) = B (0\T t (Qfah, n)... ^(xtv, t N , f N )) \0) B 

where the correspondence is valid if (ir,y) is identified with (x,it). Using asymptotic 
completeness the boundary state can be expanded in the basis of asymptotic in states as^l 

(B\ = (0|jl + 0) (2.5) 

+ I" S / Tn w J\ R2 ^ k)A m (-K, -k)A m {K, %) + ... \ 
Jo 2tt J (27i) D - 1 uj(K,k) J 

which we refer to as the cluster expansion for the boundary state (where due to translational 
invariance only bulk multi-particle states with zero total momentum can appear). 

2.2 Relation between the two channels: K 1 and K 2 in terms of R 

The one-point function of the field, due to unbroken Poincare symmetry along the bound- 
ary, only has a nontrivial dependence on x: 

B (0Mx,t,f)\0) B = Gl dry (x) 

which corresponds to a nontrivial vacuum configuration in the presence of the boundary 
condition. The leading asymptotic behaviour for x — ► — oo is given by [5] 

B (0mx,t,r)\0) B = (0m0)\0)+ge mx (2.6) 

3 The bars on top of the K coefficients indicate that the above expansion is that of the conjugate ("bra") 
boundary state. 
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where (0|$(0)|0) is the vacuum expectation value in the bulk and g is a parameter which 
is characteristic of the boundary condition (and also of the field $). We recall that \0)b is 
the ground state of the boundary system which means that there are no bulk excitations 
present and the boundary itself is in its ground state. The absence of bulk excitations 
is important for the above asymptotics to be valid; however, (12.61) also holds when the 
boundary is excited ('surface plasmons'). 

Using the property of the interpolating field $ 

(0|S(0)|A(k = fl)> = y! 

where Z is the bulk wave function renormalization constant (0 < Z < 1), and from the 
cluster expansion (12.51) one obtains the relation^ 

9 = 

On the other hand, the existence of nontrivial vacuum expectation value for the field is 
generally related to a singularity of the reflection factor at the particular kinematical point 
k = 0, uj = (i.e. k = im or equivalently •& — iir/2). In our paper [5] it was shown that 
this singularity takes the following form 

R(oj, k) ~ -^^(2n) D 5(k) (2.7) 

with g parametrizing its strength. Using the cluster property of local quantum field theory 
we proved the following relation 

- 9 [Z 

valid for general quantum field theories, which yields the expression of K 1 in terms of g: 

2 

This extends a relation previously conjectured in the case of two-dimensional integrable 
field theories (6j [T] . In the two-dimensional case, there is no parallel momentum k and the 
rapidity parametrization (12. ID takes the form 

uj = mcosh.d , n = msmhd (2-8) 

As a result, the singularity fj2.TR corresponds to a pole [TJ 

ig*/2 



$ - in/2 



4 Note that this relation remains valid if the Lagrangian field $ is replaced by any bulk interpolating 
field for the asymptotic particles and its appropriate wave function renormalization Z; in that case g also 
needs to be replaced by another constant which corresponds to the field considered. 
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Figure 3.1: The folding trick, illustrated for a generic defect scattering process 



Let us now turn to K 2 . Using the reduction formulae derived in [5] the relation to R can 
be obtained as follows: 

K 2 (k, k) = R{uj — > —in, k) 
This can be written using the rapidity parametrizations (I2.1|2.3I) as^l 



K 2 (tf, k) = R [i- + $,k 

and this relation fits very well with the pictorial representation in figure 12 . 1L In two space- 
time dimensions this is the same as the relation obtained by Ghoshal and Zamolodchiko\|§ 
in [T]. We remark that when the theory in the bulk is free and the reflection is elastic, the 
boundary state can be written in a closed forrrfl 

(B\ = (0| exp ji^A^O, 0) (2.9) 

+ l /o w r , K 2 (K,k)A in (-K,-k)A in (K,k)\ 

Jo 2-n J (27r) D ^ 1 uj(K,k) J 



3 Defects and defect operators 

Boundary conditions considered in the context of the Casimir effect generally allow trans- 
mission as well, and such boundaries are called 'defects'. A suitable generalization of the 
boundary state formalism can be obtained by a folding trick depicted in figure 13.11 which 
maps the defect into a boundary system [8]. Suppose now that a defect is located at x . 
In the crossed channel (where time flows in the x direction) it can be represented by a 
defect operator which acts from the bulk Hilbert space of the x < xo system into that 

5 Note that the rapidity arguments on the two sides of the equality are conceptually different, since they 
correspond to the kinematical variables of two different channels as defined in l|2.1|) and (|2.3|l . We can 
consider them related by analytic continuation. 

6 They also noted that the relation between the two channels can be considered as a generalization of 
the well-known crossing symmetry to quantum field theories with boundary. 

7 In 1+1 dimensions this can be extended to any integrable QFT with integrable boundary condition 
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R:>| R + :|< T + : / T": \ 
Figure 3.2: One-particle defect amplitudes 

of the x > x system^. Let us denote the operator creating the particle for the x < x 
domain as A\ while for the x > xq domain as A\. There are now four one-particle reflection 
amplitudes, shown in figure 13.21 Two of them are denoted R ± and preserve the species 
number 1,2, corresponding in the defect picture to reflections on the left and on the right 
side, respectively. The other two, T 1 * 1 are the ones changing 1 into 2 and 2 into 1, and in the 
defect picture they describe transmission from left to right and right to left, respectively. 
These can be conveniently put together into a defect matrix Q 



R+(-d,k) T-($,k) 
T+(0,jfe) R-(0,k) 



D(ti, k) = 

Using the folding map to the boundary system we obtain the defect operator [8j as^| 



° = i+ rf/(£^( fl+ (T- tf ' i >!(- tf --^< ! '-*» + < xi) 

T+(y - «, k)A\(-i), -k)A,(-0, -k) + T-(y - fl, k) t)4(«, *) + 

R~ — i?, k s jA 2 {'d, k)A 2 (—$, —k)j + terms with more than two particles 

With the same conditions as for the boundary state (trivial bulk scattering, and elasticity 
for the combined one-particle reflection/transmission amplitude) the defect operator can 
be summed up into an exponential form similar to (|2.9D . as discussed in [3]. 



4 Derivation of Casimir energy 

We now turn to the derivation of Casimir energy of a D+l dimensional scalar field x, y) 
in a domain of width L in x (for details see [HE]). Consider two defects located at a distance 
L with defect matrices D 1 and D 2 . The ground state eigenvalue of the Hamiltonian H B in 
(12 .21) can be evaluated via the partition function. Compactifying all infinite (temporal and 

8 On the two sides of the defect, the bulk theories may differ; in general, a defect can be an interface 
between very different quantum field theories (as an example one can consider the electromagnetic field in 
the presence of an interface between two drastically different physical media) . 

1 D is not necessarily unitary, since we allow for inelastic scattering processes creating and annihilating 
particles. 

9 For the sake of simplicity here we omitted possible one-particle terms corresponding to nontrivial 
vacuum configurations, but their inclusion using the folding trick is straightforward. 
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spatial) dimensions (i.e. the D extensions perpendicular to x) to circles with perimeter T 
we can evaluate the partition function in two different ways [I]: 

Z(L, T) = Tr HB e- THB = (0| D ie - LH D 2 |0> 

where H is the bulk Hamiltonian ( 12.41) in the x channel in the domain between the two 
defects and |0) is the corresponding bulk vacuum state. Inserting a complete set of bulk 
asymptotic states we obtain 

Z(L, T) = e~ LE ° (0| £>i \n) (n| D 2 |0> e - L{E "- Eo) 

n 

Normalizing the bulk ground state energy E to 0, the first few terms can be written 
explicitly as 

1 + zZ <°l D ^ % 9) % ^> 91^2 |0) e~ L ( mcff c ° sh ^+™ cff ® cosh^') 
+ 0(e" 3mL ) 

The term 1 is the contribution from the vacuum (|n) = |0)), the next term comes from 
two-particle terms in (13.11) and the higher-order corrections come from the higher multi- 
particle terms. This is a sort of cluster expansion similar to the one used in [7], valid 
for large values of the volume L. Finite volume restricts the momenta to n = and 
ki = Tjrrii, and the normalization of the creation operators becomes 

[A in ( K , k), AUk, k')} = T d u;(k, £) 

The ground state (Casimir) energy (per unit transverse area) can be extracted from the 
partition function as 



E(L) = -\im o ± 5 \ogZ(L,T) 



The result is 



-^ cM Iww^ m °" lk)x (41) 

The correction terms correspond to higher particle terms in the expansion (13. of the defect 
operator D and include the amplitudes of reflection/transmission processes involving more 
than one particle in at least one of the asymptotic states. These can be computed (together 
with the reflection factors i? ± ) e.g. using a BQFT formulation as the one presented in [3], 
but it is obvious that they are suppressed by a factor e~ mL with respect to the leading 
order term due to the presence of at least one additional particle in the corresponding term 
of the expansion of the defect operator D. Note that (I4.ip is applicable in the presence 
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of nontrivial bulk and boundary interactions: their effects at leading order are contained 
in the reflection factors R^, so as long as there is some theoretical or experimental input 
from which these can be determined the leading order contribution can be evaluated. 
In the elastic case the expansion can be summed up [4]0 

f°° dti , „ f d D ~ l k - . , 

E[L) = L^'/mk"* 1 "" (42 » 

log(l - + d,k)RZ(j _ tf,jt) e - 2 "w<*> 1 °»*'') 

Let us now calculate the ground state energy in the presence of nontrivial vacuum con- 
figuration of the field. For simplicity we suppose that the boundary is totally reflective. 
Compactifying the other directions again to circles of perimeter T with periodic boundary 
conditions we obtain 



Z(L,T) = {B a \e- LB \B p )= y £, 



-iy/m.\ - ( B <*\ n )( n \ B p) „-L'„L 
(n\n) 



-c 



The leading finite size correction to the ground state energy for large L is now given by 
one-particle terms, and the ground state energy per transverse area (at leading order in L) 
has the form [5] 

Ef{L) = -mKlKy mL + ... (4.3) 

For partially reflecting boundaries (i.e. defects) the appropriate K 1 coefficient is the one- 
particle coupling of the defects evaluated in the domain between them. If one of the K 1 -^ 
is zero then the leading correction comes from two-particle states, and is identical to (14.11) . 



5 Summary and discussion 

A very appealing property of the boundary state approach is the universality of the formulae 
04.1 p and (14.21) . In [I] we showed that the latter indeed reproduces all the results previously 
known for the planar situation, including the famous Lifshitz formula [10] (it also provides 
a way to compute new cases easily, as we demonstrated for a massive scalar field with 
Robin boundary condition). 

Another important point is that this approach formulates the Casimir effect from an 
infrared viewpoint. Standard derivations of the Casimir effect solve the microscopic field 
theory. This necessitates tackling diverse issues such as renormalization, and also the 
possibility that the infrared (long distance behaviour) may be quite different from the 
microscopic description of the theory (as is the case for example in QCD). Formula (14.11) 
expresses the effect in terms of the asymptotic particled3, and provides a long distance 
expansion for Casimir energy. 

10 We remark that the usual zero mode summation method leads to the same result, as indicated in 
Appendix A of [4]. 

1 Indeed it can be thought of as an expansion in the number of virtual particles exchanged between the 
defects. 
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Our results are consistent with the philosophy behind the more recent approach by 
Emig et. al. [llj, the origins of which can be found in the earlier papers [12l [T3J, E]. From 
this viewpoint the Casimir effect is an interaction of fluctuating surface charge densities, 
and therefore it does not logically imply the existence of (astronomically large) zero point 
energies because the bulk energy density can be trivially discarded. In the boundary state 
approach the surface is characterized by the coefficients in the cluster expansion of the 
boundary state (12.51) (or, more generally, the defect operator (13.11) ). Both approaches 
give manifestly finite results, with no ultraviolet divergences whatsoever. There are some 
differences, however. While the boundary state approach only works easily for the planar 
case, their methods can be used for general geometries. On the other hand, the approach 
of [H] is only formulated for free field theories with linear boundary conditions since it 
relies heavily on the computation of Gaussian path integrals, while in the boundary state 
approach the expansion can be written down for interacting field theories with nonlinear 
boundary conditions, in terms of their long distance scattering data. The fact that the 
path integral is Gaussian also gives Emig et al. the ability to tackle theories with zero mass 
gap, which is only possible in the boundary state approach whenever the expansion can be 
resummed into the form (14.21) . The boundary state approach, on the other hand, provides 
access to highly nontrivially interacting theories with a mass gap (a prominent example 
of which is QCD), provided the relevant scattering data are determined e.g. from lattice 
field theory (we remark that it is also highly successful in two-dimensional integrable field 
theories where exact scattering amplitudes are known). 

It is important to note that the restriction of the boundary state approach to the planar 
case comes from the fact that the high symmetry of the planar situation is exploited to 
relate the boundary states (or defect operators) to the scattering data, therefore it is not a 
restriction inherent in any theoretical principle. Finally we remark that the results (14.1114.21) 
automatically include the contribution of states localized to the defects ('surface plasmons') 
as discussed in [I]. 
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